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Abstract. The surface magnetization of the quantum XY chain in a transverse field 
is found for arbitrary nearest neighbour interactions in closed form. This allows to 
derive the bulk phase diagram in a simple way. The magnetic surface behaviour and 
the bulk correlation length are found exactly. 



We consider the ferromagnetic anisotropic XY quantum chain in a transverse field 
with hamiltonian [1] 
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The as are the Pauh matrices, /i„ is the local transverse field at site n, |7y„| < 1 the 
local anisotropy between couplings in the x and y direction. When all the r/s are 
equal to unity one recovers the Ising quantum chain [2]. With a constant anisotropy 
factor, r/„ — rj, the phase diagram and bulk critical behaviour are known exactly [3]. 
Here we focus on the surface critical behaviour of the chain with free boundary 
conditions. To diagonalize the hamiltonian, one has first to perform a Jordan- 
Wigner [4] transformation, introducing fcrmionic creation and anihilation operators 
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cj, , c„ which satisfy the usual fermionic anticommutation rules. In terms of these new 



operators the hamiltonian (1) is a quadratic form 
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where A is a symmetric matrix such that 
and B an antisymmetric matrix with elements 
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The more cfBcient way to diagonalize the quadratic form (2) is to follow the procedure 
introduced by Lieb, Schultz and Mattis [5] which is valid for a general distribution 
of the coupling constants. To proceed, we introduce a new set of diagonal fermionic 
operators e~^,e which are linear combination of the c^ and c operators: 
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where the (t>„{n) and ip„{n) are real. The hamiltonian (2) is diagonal 
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when the vectors ^ {^ ) with components i/^ {n){(f) Jn)) and the fermion excitation 
energies A satisfy the matrix relations: 



A $ 

q q 



(A-B)*^ 



(7) 



One may combine the previous matrix equations to deduce K from the eigenvalue 



problem 
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In order to solve the problem, one has to find the eigenvalues and eigenvectors of the 
excitation matrix (A — B)(A + B). All the physical quantities, such as magnetization, 
energy density or correlation functions can be expressed in terms of the eigenvectors 
$g, \l/g and corresponding excitation energies A^. 

We now concentrate on the surface magnetic behaviour. The x{y) component of 
the surface magnetization is obtained from the autocorrelation function Gs{t) = 
(a^ {0)ai (t)) in imaginary time r where a^ (t) = e'^^'^cr^'^'e^^^. In the limit 
of large r we have 

lim C^Jy^T) ^ U^(^)j^ (9) 

Using the diagonal basis of TC, we obtain 
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where |0} is the ground state and \a) = £j'"|0} is the first excited state with one fermion. 
Noticing that af — c'{ + c^ and that (j\ ~ i{c^ — c^) and making use of (5) we obtain 
m^ — 01 (1) ^iid similarly m^ — fAill)- 

Following Peschel [6] , it is now possible to obtain a closed analytical formula for the 
x{y) component of the surface magnetization. On the semi- infinite system (N -^ cxd), 
in the ordered phase the first gap E^ — Eq — A^ vanishes due to spontaneous symmetry 
breaking. In this case, equation (7) simplifies into 



(A-B)*i =0 
(A + B)$i=0 
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Noticing that changing ry into —77, m^ and m^ exchange, in the foUowing, we only 
concentrate on the x component of the surface magnetization. 

To find the eigenvetor ^^ we can rewrite the matrix equation (A + B)<i>]^ = 
under the form 



01 (n) 
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where T^ is the two by two matrix 




Iterating equation (12) one arrives at 
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The (n + l)th component of the eigenvector $-^ is then given by the expression 
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where (1|...|1) stands for the 1,1 component of the corresponding matrix. The 
normahzation of the eigenvector, ^^ <Pi{n) = 1, leads to the surface magnetization: 
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For 77 = 1 we recover Peschel's formula [6] for the Ising model in a transverse field. 
This expression is valid for any distribution of the coupling constants of the chain 
provided that the first gap vanishes, i. e., in the ordered phase(s). The transition to 
the paramagnetic phase is characterized by a divergence of the sum in (16). Since for 
a one-dimensional quantum system with short-range interactions, the surface cannot 
order by itself, the surface transition is the signal of a transition in the bulk. 

We consider in the following the homogeneous chain, 77^^ = 77, J„ = 1 and 
h^ = h \fn. In this case, all the T matrices are the same and it is easy to evaluate (16). 
The eigenvalues of T are 
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for /i^ + 77^ > 1 and complex conjugates otherwise: 
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On the line h^ + rj^ = 1 they become degenerate. First we notice that the leading 
eigenvalue gives the behaviour of 4'\{n) ^ l^^+P" which for an ordered phase implies 
that |A_|_| < 1, the first mode is then localized near the suface. From equations (17a, b) 
this corresponds to h < h^ — 1. When /i > 1 a vanishing surface magnetization is 
obtained and the bulk is in the paramagnetic phase. The line h — 1 is the Ising 
transition line since for 77 7^ it belongs to the Ising universality class. As it can be 
expected from (17) and shown hereafter, the line h^ + rj^ = 1 separates an ordinary 
ferromagnetic phase (/i^ + ry^ > 1) from an oscillatory one (h'^ + rf' < 1). This line is 
known as the disorder line of the model. Finally, the line with vanishing anisotropy 
77 = corresponds to a continuous transition (with a diverging correlation length) and 
is called the anisotropic transition line where the magnetization changes from the x- 
to the y-direction. In the ordered phase, the decay of the eigenstate ^^ gives the bulk 
correlation length which is related to the leading eigenvalue A_|_ of the T matrix by 
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so that 



e 



1 



21n|A 



(19) 



Using (17), one obtains the correlation length in both the normal and oscillating 
phases [3]. In particular, one recovers the known result that the correlation length 
exponent i^ = 1 at both the Ising and anisotropic transitions. 

Let us consider first the Ising transition line at a fixed anisotropy rj. Then, close 



enough to the transition, we have h^ 
easy to show that 



77 > 1. The eigenvalues A_(_ are real and it is 
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where AA = A^ — A_. This expression reduces for large n to A"^^/AA. Using 
equation (16) we obtain the surface magnetization near the Ising transition 
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Thus the surface critical exponent is /3f = 1/2 on the Ising transition line. 

Close enough to the anisotropic transition line (77 — 0), we are in the oscillatory 
phase. The complex conjugate eigenvalues A_(_ of equation (17b) can be parametrized 
as 
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with 
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Using (14) a straightforward calculation shows that 



0iH = (-!)" 



l + ?7 



Vl - ??2 - h^ 



A" -A" 



M^) 



4)",".i±4sh 



(-l)>"/n^^-("^)^^(l) 



(24) 



We sec from this expression the oscillatory nature of the phase with a wavelength 
given by ^^^[3]. Close to the transition the normalization of the eigenvector leads to 
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Finally one can also consider the behaviour on the disorder line h"^ + if' = 1 
separating the oscillatory phase from the normal one. On this line, the eigenvector 
reads 
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with p defined in (23). Thus we obtain 
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Near the special point (77 = 0, /i = 1) we have simply rrvl ~ 2rf''^ with a surface 
magnetic exponent equal to 3/2. This special value is due to the merging of the 
two relevant fields of the model. We may notice that this expression can be deduced 
from (21) or (25) using the constraint h? + rf' = 1. 

To summarize, we have determined the surface magnetic behaviour of the XY 
quantum chain using a very simple analysis. We have shown how to deduce the 
bulk phase diagram and the correlation length from the surface behaviour. This 
procedure may be of interest for other one-dimensionnal quantum systems where the 
surface behaviour is more tractable than the bulk one. Finally one may notice that 
expression (16), cut off at some size L, may be used for finite-size scaling calculations 
with arbitrary couplings [7]. 
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